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Cosmological phase transitions are predicted by Particle Physics models, and have 
a variety of important cosmological consequences, which depend strongly on the dy- 
namics of the transition. In this work we investigate in detail the general features 
of the development of a first-order phase transition. We find thermodynamical con- 
<^ [ straints on some quantities that determine the dynamics, namely, the latent heat, the 

radiation energy density and the false- vacuum energy density. Using a simple model 
with a Higgs field, we study numerically the amount and duration of supercooling 
qj . and the subsequent reheating and phase coexistence. We analyze the dependence 

\ of the dynamics on the different parameters of the model, namely, the energy scale, 

0^ \ the number of degrees of freedom and the couplings of the scalar field with bosons 

and fermions. We also inspect the implications for the cosmological outcomes of the 
phase transition. 

43 , . I. INTRODUCTION 

Particle Physics models predict the occurrence of several phase transitions in the early 
Universe, such as e.g., the electroweak phase transition or the quark-hadron phase transition. 
Phase transitions in the early Universe may leave observable vestiges, such as topological 
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rnase transitions in trie early universe may leave observable vestiges, sucn as topological 
defects pL magnetic fields jij, the baryon asymmetry of the Universe jjjj], baryon inhomo- 
geneities |3, H|, gravitational waves 0] or black holes Q. The effects of some of these relics 
can constrain the model, as in the case, e.g., of monopoles and domain walls. Late time 
cosmological phase transitions have also been proposed to act as seeds of the large-scale 
structure formation and as an explanation of the dark energy problem [§, 0, 10, The 
outcome of a phase transition depends, both quantitatively and qualitatively, on several 
^ ■ aspects of the dynamics, for instance, the nucleation rate, the velocity of bubble expansion, 
and the temperature variation during the development of the transition. 

In general, the evolution of a first-order phase transition can be divided in three stages, 
namely, supercooling, reheating and phase coexistence. At T = T c the free energy has two 
degenerate minima separated by a barrier. Hence, the bubble nucleation rate T vanishes. At 
T < T c , "critical" bubbles nucleate. These are bubbles of the true vacuum which are large 
enough that their volume energy dominates over their surface tension, so they can expand. 
Assuming that the standard picture of bubble nucleation applies, bubbles of the supercooled 
phase will nucleate in a homogeneous background of true vacuum. The number of bubbles 
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will not be appreciable until a lower temperature T/v, which can be estimated as follows. 
The age of the Universe is t ~ H~ l , and a causal volume is Vjj ~ t 3 . Then, if at least one 
bubble is to be created in a time ~ t in a volume ~ Vh, we must require that Tt 4 > 1. Thus, 
the temperature T/v is roughly determined by the condition T ~ iT 4 . 

In fact, this picture may not work and the supercooling stage may be shorter (or not occur 
at all). For instance, the presence of impurities (such as e.g. topological or non-topological 
solitons) could trigger bubble nucleation [12]. Also, if the phase transition is weakly first- 
order, i.e., if the barrier of the free energy is sufficiently small, thermal fluctuations called 
subcritical bubbles may dominate [l3[. In this case, there may be a two-phase emulsion 
already at T = T c . Then, sub-critical bubbles may percolate and true-vacuum domains may 
begin to grow at a temperature T > T N . 

Initially, bubbles of true vacuum grow with a velocity which is governed by the pressure 
difference across their walls and by the viscosity of the hot plasma or relativistic gas sur- 
rounding them. As bubbles expand, latent heat is liberated and reheats the system back 
to a temperature T r . As a consequence, the expansion of bubbles slows down, since the 
pressure difference decreases as T approaches T c . If the latent heat L is negligible, there will 
be no temperature variation. One expects that reheating will be important if L provides 
the energy density difference needed to increase the temperature of radiation from T/v back 
to T r ps T c , i.e., when L ~ Spn ~ T c 4 — T^. If L is much larger than 8p Rl the temperature T r 
will be very close to T c . When this happens, a stage of "slow growth" or "phase coexistence" 
follows. Indeed, since T cannot increase beyond T c , bubbles will grow only at the rate at 
which the expansion of the Universe takes away the injected energy. The temperature will 
thus remain nearly constant until every region of space has been converted to the stable 
phase. 

Although the above picture is quite general, the details of the dynamics depend on the 
specific model. A complete analysis involves, even in the simplest cases, solving a set of 
integro-differential equations for the nucleation and expansion of bubbles, which takes into 
account the reheating of the thermal bath. Therefore, it is useful to find general char- 
acteristics, which will permit to obtain some conclusions before embarking on the task of 



computing the development of a given phase transition. In Ref . [14j , an analytical approach 
was performed, which allowed to obtain some general conclusions on the evolution. However, 
due to the involved dynamics of reheating, the analytical study requires some rough approx- 
imations, particularly for the nucleation rate. A numerical investigation is thus necessary 
in order to have a better understanding of the dynamics of first-order phase transitions and 
their cosmological consequences. 

In this work, we shall perform a detailed study of the general dynamics of phase transi- 
tions. We shall be interested in first-order phase transitions occurring either in the radiation 
dominated epoch, or in a sector composed of radiation. In particular, we shall examine ther- 
modynamic constraints which apply to any first-order phase transition. As we shall see, this 
allows to discuss on the possible effects of a model without making numerical calculations. 
We shall also make a numerical investigation of the dynamics. For that purpose, we shall 
use a simple model for the free energy, which allows to consider different kinds of phase 
transitions, both weak and strong. The model also provides an approximation for realistic 
theories (e.g., different extensions of the Standard Model). We shall discuss the implications 
of our results for the cosmological outcomes of the phase transition. 

The article is organized as follows. In the next section we discuss some general properties 
of phase transition dynamics, and we study model-independent relations between thermody- 



namical parameters. Then, in section II I II we consider a simple model, consisting of a scalar 
(Higgs) field, which has Yukawa couplings to different species of bosons and fermions. We 
write down the one-loop finite-temperature effective potential for this model, and discuss 
the different kinds of phase transitions the model can present. In section IIVI we consider 
the equations for the evolution of the phase transition, and we compute them numerically. 
We are particularly concerned with the amount and duration of supercooling, and with the 
extent of the phase coexistence stage. 

We apply the results of this investigation in section |V] where we analyze some of the 
possible cosmological outcomes of a phase transition to illustrate the effect of the dynamics. 
We consider the formation of baryon inhomogeneities in the electroweak phase transition, 
the creation of topological defects, and the generation of magnetic fields. We also discuss on 
different proposals of late-time phase transitions as solutions to the dark-energy problem. 
We show that thermodynamical constraints rule out some of these models. Our conclusions 
are summarized in section I VI I Some technical details of the calculation are left to the 
appendix. 

II. PHASE TRANSITION AND THERMODYNAMIC PARAMETERS 

We can use thermodynamic considerations to obtain some general information on the 
amounts of supercooling and reheating and on the duration of the phase transition, without 
specifying the form of the free energy. 

A. Supercooling and phase coexistence 

Consider a system which undergoes a phase transition at a temperature T c . The high- 
temperature phase consists only of radiation and false vacuum energy, so the energy density 
is of the form 

P+ = PA + PR, (1) 

where p\ is a constant and pn = g*7r 2 T 4 /30, where g* is the number of relativistic degrees 
of freedom (d.o.f.). At the critical temperature T c the two phases have the same free energy 
density, but different energy density. The discontinuity Ap(T c ) = p + {T c ) — p_(T c ), with p^ 
the energy density of the low-temperature phase, gives the latent heat 

L = Ap(T c )=T c As(T c ) } (2) 

with As = s + — s_ the entropy density difference. This entropy is liberated as regions which 
are in the high-T phase convert to the low-T one. 

Since entropy is conserved in the adiabatic expansion of the Universe, the entropy density 
of the system can be written as 

s = s + (T c )(ai/a) 3 , (3) 

where a is the scale factor, and <2j is its value at the beginning of the transition, i.e., at 
T = T c . During the phase transition, s is given by 

s = s + (T)-As(T)f, (4) 



where / is the fraction of volume occupied by bubbles of low-T phase. 



If there is little supercooling (e.g., if the phase transition is weakly first-order, or if bubble 
nucleation is triggered by impurities), the temperature T N at which bubbles form and start 
to grow will be very close to T c . In this case, a small L can take the system back to T c . 
Then, a good approximation is to consider that the phase transition develops entirely at 
T = T c , with equilibrium of phases [RlS]. Thus, the fraction of volume is easily obtained 



from Eqs. (3) and @). The result is [3 

/ - s+(7; 



a . 



(5) 



As(T c 

The phase transition completes when / = 1, so its duration is determined by the condition 

(oi/aff = l-As(T e )/s + (T c ), (6) 

where a/ is the scale factor at the end of the phase transition. 

In general, though, bubble nucleation does not begin as soon as T reaches T c . The 
temperature decreases until the nucleation rate becomes comparable to the expansion rate. 
During supercooling, the entropy of the system is that of radiation, s+(T) = sr(T), with 

s R (T) = = ^T 3 , (7) 
3 T 45 K ' 

so, from Eq. ([3]) we have T = T c a,i/a. When the number of bubbles becomes noticeable, the 
released entropy begins to reheat the system. The minimum temperature T m delimits the 
end of supercooling. It is reached at a value a m of the scale factor given by T m « T c ai/a m . 
One expects that for L > 5p R = Pr(T c ) — pit(T m ), the temperature will go back to T « T c 
and a period of phase coexistence will begin. We will now show that the condition for phase 
coexistence to occur is in fact 

As(T c ) > 6s R , (8) 

where 5sr = sr(T c ) — s#(T TO ). In terms of energy, we have T c As(T c ) = L and 5sr = 
(4/3)5 (pr/T), so the above condition becomes L > (4/3) SpR. 

Assuming that a phase coexistence stage at T r w T c is reached, we can go back to Eqs. 
(J3J) and (02), which lead again to the result (jSJ) for the total change of scale a//aj, even though 
this time the temperature was not constant from the beginning. Therefore, the final value 
of the scale factor a/ is not affected by the previous supercooling and reheating stages. This 
will only be possible, however, if a m < df, since the supercooling stage cannot be longer 
than the total duration of the phase transition. During supercooling, s = s+, so s + (T m ) 
is given by Eq. © with a = a m . Comparing with Eq. 0, the condition a m < cif gives 
s + (T m ) > s + (T c ) — As(T c ). Since s + = sr, Eq. (jSJ) follows. 

The value of L can be easily calculated for any model, since it is derived directly from 
the free energy. In contrast, calculating 5p + entails the evaluation of the nucleation rate T, 
which must be calculated numerically, and then solving the equations for the evolution of 
the phase transition in order to determine T m . We will perform such calculation in section 
HVl Provided that condition (jSJ) is fulfilled, the value of a/ will be independent of the amount 
of supercooling, and given by Eq. (jSJ). We can write equivalently 

(ai/a f ) 3 = l-SL/Ap R . (9) 

How long will the phase transition go on, depends on how large L is. Since the entropy 
difference is bounded by As(T c ) < s + (T c ), the latent heat has a maximum value L max = 



T c s + = 4p fi /3. We see that aj — > oo in this limit. This is because s_ = 0, so all the 
entropy must be extracted from the system in order to complete the phase transition, and 
this requires an infinite amount of work. 

The duration At of the phase transition is related to the expansion factor af/cii through 
the expansion rate H. Consequently, it depends on the different kinds of energy (e.g., matter, 
vacuum, radiation) that make up the total energy density p. If our system is uncoupled from 
other sectors (as in the case of late-time phase transitions), then it is not straightforward to 
calculate At. In the early Universe, instead, we can assume that all particle species are in 
equilibrium with each other and constitute a single system which is dominated by radiation. 
Then, for the period of phase coexistence at T = T c , the equation of state is especially 
simple, since temperature and pressure are constant. The energy density is given by 
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where 



Pc = Pr(T c )/3 - p A 
is the pressure at T = T c . Consequently, the Friedmann equation 1 
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where G is Newton's constant, can be solved analytically 15l . [16] . We have 
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where u = \J QnGp c and 5 = arcsin ^p c /T c s + . 
From Eqs. (EI) and ( {TBI we obtain 



At 
t 




arcsin 




y/l ~ As/S 



(14) 



where t = (2Hi) 1 ~ ti. Notice that At ft depends only on the two parameters r = L/ pn 
(equivalently, As/s + ) and R = p\/ Pr. We remark that, as long as a temperature T r w T c 
is reached after reheating, At gives the total duration of the phase transition, i.e., the time 
elapsed from the beginning of supercooling at t — ti until the end of phase coexistence at 
t = tj. As we have seen, the condition for the validity of Eq. (|14p is that supercooling ends 
before this time. Otherwise, At will be given essentially by the duration of supercooling, 
since the subsequent reheating stage will be short. In that case, Eq. f|T4l) gives just a lower 
bound for the duration of the phase transition. 



B. Constraints on thermodynamic parameters 

At the critical temperature, one expects that the energy density of radiation is at least 
of the order of that of the false vacuum, since radiation must provide the entropy necessary 
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to make the minima of the free energy degenerate. Notice that the exact relation between 
Pr(T c ) and p\ can be determinant for the dynamics of phase coexistence. Indeed, for p R /3 < 
Pa, the pressure p c is negative and the sine in Eq. (1131) becomes a hyperbolic sine, which 
indicates that the expansion of the Universe is accelerated. This happens because the energy 
density (TTOT) includes a constant term p A ff = — p c , which represents an effective cosmological 
constant |9|, ll7|. In this case, p A ff > 0. Moreover, if pa ~ pr(T c ), the false vacuum energy 
may become important before the phase transition, i.e., at T > T c . 

On the other hand, if pa < P_r/3, we have p A ff = — p c < 0. Then, according to Eq. ffT3]) 
the Universe will collapse after a time t c ~ 1/w, unless phase coexistence ends before this 
time, so that this equation is no longer valid. Notice that phase coexistence may be long 
if L pa 4p#/3. The collapse occurs because the energy density ffTUj) and, consequently, the 
expansion rate ( TT2l) vanish for a finite value of a/cij. Nevertheless, the quantities L, pr, and 
Pa are constrained by thermodynamical relations, and we will show that none of the above 
situations can arise, i.e., phase coexistence will not cause either accelerated expansion nor 
collapse of the Universe. 

The pressure of the relativistic system is given by p = —J 7 , where T is the free energy 
density. Hence, at T = T c we have p c = —T+ (T c ) = — JF_ (T c ). The free energy density 
depends only on temperature, dT = —sdT — (p + JF) dV/V = —sdT. Since s > 0, T (T) 
must be a monotonically decreasing function. Therefore we have in particular JF_ (T) < 
T- (T = 0) for any T > 0. But at T = the free energy matches the energy. Hence, assuming 
that the energy density vanishes in the true vacuum, we have JF_ (T = 0) = p-(T = 0) = 0. 
Then, JF_ (T c ) < and p c > 0, so the condition for accelerated expansion is never fulfilled. 
Moreover, the condition 

PA < Pr(T c )/3 (15) 

implies that false vacuum energy never dominates, unless the system departs from thermal 
equilibrium (for instance, pa may become dominating in the course of supercooling). 

Now, since p c > 0, the Universe will not collapse only if phase coexistence ends before p 
vanishes. According to Eq. (fTUj) . this is true if (ay/cij) 3 < T c s + /p c . Using Eqs. ©, © and 
(|TT|) the condition becomes 

L < p + = p A + Pr{T c ). (16) 

But this is always fulfilled, since L = p + — p_, and P-(T) > at T > [because dp/dT = 
Tds/dT> and p_(T = 0) = 0]. 

The inequalities above become equalities only for JF_(T) = p~{T) = 0, i.e., at T = 0. So, 
both limiting values pa = Pr/3 and L = pa + Pr are attained only if T c = 0. In this limit pa 
and pr vanish, but still L/p R 4/3. Hence, Eq. implies that af — * oo. Thus, for a phase 
transition with T c pa we will have a long phase-coexistence stage. For a given model with 
a fixed energy scale v , small T c means T c v, i.e., the metastable minimum and the barrier 
must persist at T v. At such low temperatures, the free energy coincides approximately 
with the zero-temperature potential, and the minimum <p c tends to the zero-temperature 
value v. This corresponds to a very strongly first-order phase transition, with C /T C ^> 1. 
In this case one expects that the nucleation rate will be suppressed and the supercooling 
stage will be long too. However, it is not straightforward to compare the duration At s of 
supercooling to that of phase coexistence, since the latter depends significantly on the total 
number of d.o.f. g*, while At s depends essentially on the bubble nucleation rate T. In 
section [IV] we will see that, depending on the model, we can have either At s At (i.e., 
little supercooling) or At s pa At (i.e., short phase coexistence). 



In a specific model, the parameters p\, L, and pr can be derived from the free energy. The 
constraints (11511161) . i.e., R < 1/3 and r < R + 1, should then be automatically fulfilled 2 . In 
general, pa and L will be even more constrained. For instance, the radiation density pr may 
contain a component pi from particles which are in thermal equilibrium with the system, but 
are not directly coupled to the order parameter, and therefore do not contribute to L and 
Pa (e.g., "light" particles which do not acquire masses through the Higgs mechanism). The 
inequalities above hold for the radiation of the system alone, i.e., Pr — pi, so the constraints 
become pa < (pr — pi) /3 and L < p\ + p R — p\ . If pi = ^7r 2 T 4 /30, we have R < x/3 and 
r < R + x, where x = 1 — g\j g*. 

Fig. [TJ shows the region in the (R, r)-plane allowed by thermodynamics, and inside that, 
the contours of constant time At. On the right axis we have indicated some values of 
af/a.i (which depend only on r). The points correspond to some of the phase transitions 
considered in the next section. We have plotted two sets of curves, corresponding to g\ = 
and gi/g* ~ 0.44. The dashed line delimits the allowed region for the latter case. As the 
phase transition becomes stronger, the latent heat increases. However, the limit L = pa + pr 
is reached for T c — > 0, together with the limit pa = P_r/3. That is why all the curves approach 
the upper-right corner of the allowed region. 
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FIG. 1: Contours of constant time in the allowed region of the plane (p\/ Pr,L/ pr). From bottom 
to top, the curves correspond to At/t = 0.2, 0.5, 1, 2, 3 and 5. The points correspond to varying /i^ 
in the model of section ITTT1 for hf = 0.7 (blue squares), hf = hb (black triangles), and hf = hf, with 
Pb 7^ (red circles). The three curves on the right correspond to g\ = 0, and those on the left to 
gijg* w 0.44. 

The analytic approximation given by Eq. f|T4l for the total duration of the phase tran- 
sition is valid only if condition ([H|) is satisfied. Furthermore, we cannot describe, within 
this approach, the transition between supercooling and phase coexistence, i.e., the reheating 
stage. A complete description of phase transition dynamics involves the computation of the 
nucleation rate. This requires specifying a model for the free energy. 



2 Notice that some approximations for the free energy may allow values that fall outside this region (see 
e.g. the discussion on dark-energy models in section fV|) . 



III. THE FREE ENERGY 



We will consider a theory described by a scalar field with tree-level potential 



M0) = -^V+^> 4 , 



(17) 



which has a maximum at = and a minimum at = v. The one-loop effective potential 
is of the form 

V( ( p) = V ( ( p) + V 1 ( ( p)+p A , (18) 

where V\ (0) is the one-loop zero-temperature correction, and we have added a constant 
Pa so that the energy density vanishes in the true vacuum. Imposing the renormalization 
conditions that the minimum of the potential and the mass of do not change with respect 
to their tree-level values [13], the one- loop correction is given by 
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where g>j is the number of d.o.f. of each particle species, m 8 (0) is the 0-dependent mass, 
and the upper and lower signs correspond to bosons and fermions, respectively. 

The free energy density results from adding finite-temperature corrections to the effective 
potential, 

^(0,T) = K(0)+^(0,T), (20) 
where the one-loop contribution is 
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and J_, I + stand for the contributions from bosons and fermions, respectively, 
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For simplicity, we will consider in general masses of the form rrii (0) = hi<f>, where hi is the 
Yukawa coupling. Thus, the free energy takes the form 



^(0,T) = V (<j>) + J2 
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where the last term accounts for the contribution of species with hi — 0, so g\ is the effective 
number of d.o.f. of relativistic particles. The constant pa is obtained by imposing that 
V(v) = 0, so 

32vr 2 / 4 ' 
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Notice that p\ gives the energy density of the false vacuum, p\ = V (0) . 



At high temperature the free energy (1231) has a single minimum at <fi = 0. As the 
temperature decreases, a non-zero local minimum <p m (T) develops. Therefore, the free 
energy in the high- and low-temperature phases is given by J- + (T) = T (0, T) and T- (T) = 
J- {4>m (T) , T), respectively. In the phase with = 0, all particles are massless and 

F + = -g^ 2 T 4 /90 + p A , (25) 

where g* = J2db + (7/8) 9} ls the effective number of d.o.f. {b stands for bosons and / 
for fermions). Thus we have radiation and false vacuum. At the critical temperature T c , the 
two minima = and m (T c ) = C have the same free energy. Below this temperature, 
m (T) becomes the global minimum. In general, as temperature decreases further the 
barrier between minima disappears and the minimum at = becomes a maximum. This 
happens at a temperature T given by 

T 2 _ A + E Tgihj/I67i 2 2 
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Finally, at zero temperature we have JF(0, 0) = V((p), so m (O) = v. Notice, however, 
that the zero-temperature boson contribution may turn the maximum at = of the tree- 
level potential Vo(0) into a minimum of V(4>). In this case, there will be two minima still at 
T = 0. Indeed, for J^gbhf > ^ gfhj + 167i 2 X, the r.h.s. of Eq. ( f26l) becomes negative, which 
means that the barrier never disappears. Furthermore, for strongly coupled bosons the origin 
can become the stable zero-temperature minimum. Indeed, for J2gbhf > J2dfh} + 327r 2 A 
the vacuum energy density fflM becomes negative. In that case, the origin is stable at all 
temperatures, and there is no phase transition. 

The energy density can be derived from the free energy by means of the relations p = Ts + 
T and s = —dJ r /dT. Thus, from Eq. (|25l) we obtain p + = pa + pn, and p_ = —TT'_ + 
At T = T c , JF + = JF_, so the latent heat is L = — T c AjF'. Taking into account that 
dT / 8^=^ = 0, we find 
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The functions I± (x) are negative and monotonically increasing, so we see that the one-loop 
effective potential satisfies the thermodynamical bound L < —2g{T^l T (0) /ir 2 = 4p^/3. 
Furthermore, I + and J_ fall exponentially for large x. Therefore, L approaches the limit 
L/ p R — > 4/3 for hi<p c /T c — > oo. 

For our purposes it will be sufficient to consider only four particle species, namely, two 
bosons and two fermions. In this way we can have weakly coupled fermions and bosons 
with Yukawa couplings hfi and h bh and d.o.f. gji and g^, respectively. These particles 
will be relatively light in the low-temperature phase. We will consider also g b bosons and gf 
fermions with variable couplings hb and hf, respectively. The values of the Yukawa couplings 
are constrained by perturbativity of the theory, which sets a generic upper bound hi < 3.5 
Till. In addition, we include gi light d.o.f., for which we assume hi = 0. This model allows 



us to explore several kinds of phase transitions. 

For instance, choosing v = 2A6GeV we have a phase transition at the electroweak scale. 
We obtain a good approximation for the free energy of the Standard Model (SM) if we 
consider gfi = 12 fermion d.o.f. with hfi ~ 0.7 (corresponding to the top), and gu = 6 boson 
d.o.f. with hu ~ 0.35 (corresponding to the transverse gauge vectors W and Z). The rest 



of the SM d.o.f. have h <C 1, so their contribution to the ^-dependent part of the effective 
potential is negligible. They only contribute to g\ in Eq. (1231 . To make the electroweak 
phase transition strongly first-order we need to add some extra particles to the SM. For our 
purposes, we don't need to refer to any specific extension of the model. We choose A ~ 0.12, 
which corresponds to a Higgs mass m# = 120GeV, and we consider for the time being 
adding equal numbers of bosons and fermions, with gb = 9f = 10 and hb = hf = 0.7, which 
give a value <fr c /T c ~ 1.3 for the minimum of the potential at the critical temperature (see 
Fig. E]). 



0.03 
" 0.02 

K 0.01 



„ o.oo 

-0.01 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

4>/T 

FIG. 2: The free energy around T = T c . 

It is well known that heavy bosons enhance the strength of the phase transition. We 
see in the left panel of Fig. [3] that the minimum <p c , as well as the height of the barrier, 
increase if we increase the value of hb- Besides, the critical temperature decreases. Indeed, 
for fixed hf, according to Eq. (1261) the temperature T vanishes for a value hb = hbi given by 
dbhfi — ^-^ 2 ^+9fh A f+gfih A f l —g b ihl l . At this point, a barrier appears in the zero-temperature 
effective potential, and = becomes a local minimum of V(<j)). If hb is increased further, 
the zero-temperature barrier increases as the energy p\ of the origin decreases. According to 
Eq. (I24j) . the two zero-temperature minima become degenerate for a value hb = h b2 given by 
9bh\ 2 — 327r 2 A + gfhj + gfih^ — guh^. For this value of hb the critical temperature vanishes 
and 0c = v. Beyond the value hb2 there is no phase transition. 

If we keep hf = hb as we increase hb, the behavior is quite different, since the fermions 
compensate the effect of the bosons. As we can see in Fig. [3] (right panel), the height of the 
barrier increases more slowly with hb, and the value of C /T does not change significantly. 
With this symmetric choice of parameters, the false vacuum energy density (I2~4"|) does not 
depend on hb and hf. Thus, the origin will never be the stable minimum at T = 0, and T c 
will never vanish. According to Eq. (I26I) . the temperature T does not vanish either, but 
it decreases as 1/hb- In the next section we will analyze the effect of these two opposite 
variations. 

For large couplings, the first-order phase transition becomes stronger and the latent heat 
(1271) increases. The maximum value L/p^ = 4/3 will be achieved when all the couplings 
hi are large. If gi ^ 0, this maximum becomes 4x/3, with x = (1 — gi/g*)- Consider the 
case gu — gp — and gb = 9f = 10. For hf fixed, the maximum is reached at h b = h b2 , 
i.e., when T c = 0. For the case hf = 0.7 we obtain the points in the (R, r)-plane that are 
shown in blue squares in Fig. [TJ For hf = h b , on the contrary, there is no such limit on hf,. 





In this case (black triangles in Fig. [Tj), as the coupling is increased the points accumulate 
near the point (x/3, 4x/3), which is the corner of the thermodynamically allowed region. In 
particular, for gi = 0, we see that L/p R can be very close to the maximum 4/3, even though 
in this case T c does not vanish. 

It is interesting to consider the case hf > hf,. Notice, however, that strong fermion 
couplings hf may destabilize the zero-temperature potential, since they introduce negative 
quartic terms in V (0). To stabilize the potential in the case of a strongly coupled fermion, we 
can add a heavy boson with the same coupling hf and d.o.f. and a mass ml((p) = fif+hlcp 2 
ll. If \i h is large enough, this boson will be decoupled from dynamics at T ~ v. The 



maximum value of fib consistent with stability is obtained by requiring the quartic term to 
be positive for <ft 3> v. It is given by 
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(28) 



For a weakly-coupled fermion, /if, is much larger than v and the stabilizing boson is com- 
pletely decoupled. On the contrary, for large hf, m b approaches rrif and we recover the 
previous case. We have plotted in Fig. [T]the points of the (R, r)-plane 3 corresponding to a 
variation of hf (red circles). For small values of hf we have only the fermion contribution, 
and the phase transition is weakly first-order. In fact, there is a minimum value of hf for 
which the phase transition becomes second-order. At this point, the latent heat vanishes for 
a finite value of p\- In contrast, for large hf we have, as in the previous cases, a strongly 
first-order phase transition. 

As we see in Fig. HJ in all the cases the total duration of the phase transition becomes 
significant for large hi. However, the durations of supercooling and phase coexistence can 
be extremely different in each case. 



3 The definitions of p\ and pn change slightly in this case. 



IV. THE PHASE TRANSITION 



A. Phase transition dynamics 



The nucleation and growth of bubbles in a first order phase transition has been extensively 



studied (see e.g. 0, [l_8|, |20j, l2l|, |22j, [2j]). According to the conventional picture of bubble 
nucleation, at T > T c the field takes the value = throughout space. At T < T c , bubbles 
of the stable phase (i.e., with the value <ft — 4>m inside) nucleate. We remark that in a weakly 
first-order phase transition this picture may not work |13j. A quantitative determination of 
the importance of subcritical bubbles requires in general numerical calculations and is out 
of the scope of the present investigation. For instance, lattice calculations for the case of 
the minimal standard model (with unrealistically small values of the Higgs mass) show that 
subcritical bubbles may play a significant role at the onset of a weakly first-order electroweak 



phase transition [25| . Thus, our results for the amount of supercooling become unreliable in 
the limit of very small values of the coupling hb- 

The thermal tunneling probability for bubble nucleation per unit volume per unit time 



is 



27 



r ~ A (T) e~ S3/T . (29) 

The prefactor involves a determinant associated with the quantum fluctuations around the 
instanton solution. In general it cannot be evaluated analytically. However, the nucleation 
rate is dominated by the exponential in (I2"9"j) . so we will use the rough estimation A (T) ~ 
The exponent in Eq. ( l29l) is the three-dimensional instanton action 



= 4vr 



r dr 



(30) 



where AjF (0, T) = T ((f), T) — T (0, T). The configuration (r) of the nucleated bubble may 
be obtained by extremizing this action. It obeys the equation 

dT 



d 2 (p 2 d(p 
dr 2 r dr 



(31) 



Hence, S 3 coincides with the free energy that is needed to form a bubble in unstable equilib- 
rium between expansion and contraction. At the critical temperature the bubble has infinite 
radius, so S3 = 00 and T = 0. In contrast, at T = T the radius vanishes, so 5*3 = and 
T ~ T£, which is an extremely large rate in comparison to H 4 ~ (T 2 /Mp) 4 . Therefore, the 
number of bubbles will become appreciable at a temperature which is rather closer to T c 
than to To. Thus, in order to have supercooling at T < T c , the temperature To must not 
exist, so that the barrier between minima persists at T = 0. 

After a bubble is formed, it grows due to the pressure difference at its surface. There 
is a negligibly short acceleration stage until the wall reaches a terminal velocity due to the 
viscosity of the plasma (see, e.g., [28(). The velocity v w is determined by the equilibrium 
between the pressure difference p_ — p + = T + — T- = AjF(T) and the force per unit area 
due to friction with the surrounding particles, /friction = V v w- Thus, 



v w {T) = AF{T)/r). 



(32) 



The friction coefficient can be written as rj = rjTa, where rj is a dimensionless damping 



coefficient that depends on the viscosity of the medium, and a 
wall tension (for a review and a discussion see [3]). 



J (d(j)/dr) dr is the bubble 



We will assume that the system remains close to equilibrium, which is correct if v w is 
small enough. If the wall velocity is lower than the speed of sound in the relativistic plasma, 
c s = v/T/3, the wall propagates as a deflagration front. This means that a shock front 
precedes the wall, with a velocity v a h > c s . For v w <C c s , the latent heat is transmitted 
away from the wall and quickly distributed throughout space. We can take into account 
this effect by considering a homogeneous reheating of the plasma durin g th e expansion of 
bubbles 0, 29]. (For detailed treatments of hydrodynamics see, e.g., 20.l2ll|). 



The radius of a bubble that nucleates at time if and expands until time t is 

*(*-,«) = * cn ;fg+jfMn^*-. (33) 

The scale factor a takes into account the fact that the radius of a bubble increases due to 
the expansion of the Universe. The initial radius Ro can be calculated by solving Eq. ( I3TT) 
for the bubble profile <j>{r). It is roughly ~ T _1 . Hence, Rq can be neglected, since the 
second term in Eq. ( |33l) . which is determined by the dynamics, depends on the time scale 
St ~ H- 1 ~ Mp/T 2 . 

The fraction of volume occupied by bubbles is given by 

/ (t) = 1 - exp j- jf 3 r (T) ^R (f , t) 3 dt'\ . (34) 

The integral in the exponent gives the total volume of bubbles (in a unit volume) at time t, 
ignoring; overlapping. The complete expression (|34|) takes into account bubble overlapping 
30| . The factors of a take into account that the number density of nucleated bubbles 
decreases due to the expansion of the Universe. 

To integrate Eq. ([Ml) , we still need two equations in order to relate the variables T, a 
and t. Eqs ([5]) and (13) give the relation Jjl 



T ~ 27TV45 f «3 ' (35) 

where the first term, which is proportional to the released entropy As(T)f, accounts for 
reheating, and the second term accounts for the cooling of the Universe due to the adiabatic 
expansion. Finally, the Friedmann equation ( 1121) gives the relation 



Ida / 8nG 

aTt = y— p > (36) 

with p = p+ (T) - Ap (T) /, where p+ = p A + #*7r 2 T 4 /30, and Ap = -TAP + AT. 

The functions AjF(T) and AjF'(T) are easily obtained by numerically finding the mini- 
mum m (T). The nucleation rate T(T) can be calculated by solving numerically Eq. ( I3TT) 
for the bubble profile, then integrating Eq. (|30|) for the bounce action, and using the result 
in Eq. (1291) . We solved Eq. (131 p iteratively by the overshoot-undershoot method 4 . The 
thermal integrals (!2"2"j) for the finite-temperature effective potential can be computed numer- 
ically. However, we find that the computation time is lowered significantly by using instead 
low-x and high-x expansions for I±(x) (see the appendix). 



4 We have checked our program by comparing with the results of Rcf. 22 1 for the bounce and Ref. [23| for 
the evolution of the phase transition. 



B. Numerical results 



We begin by considering a phase transition at the electroweak scale, with the free energy 
plotted in Fig. El The development of the phase transition depends on the specific heat of 
the thermal bath, i.e., on the total number of d.o.f. We can take into account the light d.o.f. 
of the SM by setting g\ pa 90 in the last term of Eq. ( 1231) . The friction coefficient fj depends 
on the model and its computation is not straightforward. For the time being, let us assume 
fj ~ 1. The solid curve in Fig. H] shows the temperature variation during the phase transition 
for this model. We observe a considerable reheating, which indicates that the latent heat is 
comparable to the energy density SpR needed to take the radiation back to T = T c . However, 
a phase coexistence stage is not achieved, which reveals that L < (4/3)5p#. Notice that the 
large number of d.o.f (g* ~ 100) makes the energy density of radiation much larger than the 
latent heat. For lower values of gi, the thermal bath has a smaller specific heat and is more 
easily reheated. This can be seen in the dashed and dashed-dotted lines in Fig. HI 
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FIG. 4: The temperature variation for the potential of Fig. [2j 

We note that supercooling finishes at a temperature which is quite closer to T c than 
to T . As mentioned before, this fact is quite general, as it is due to the extremely rapid 
variation of the nucleation rate, which becomes T ~ T 4 3> H A at T = To. Notice also that 
the different curves in Fig. H] coincide during supercooling. This is because in this stage 
the relation between the dimensionless variables T/T c and r = (t — ti)/ (2H i )~ 1 is almost 
independent of any parameter of the model. Indeed, during supercooling T/T c = at/ a, and 
the dependence of the scale factor on time is given by da/a = Hdt = (1/2)(H / H^dr. Since 
a pa ai and H pa Hi for T pa T c , we have d(T/T c ) pa —{\/2)dr as long as T does not depart 
significantly from T c . 

We can check the approximation ( fl4l) for the total duration of the phase transition. The 
relevant parameters r = L/pR and R = Pa/pr are different for each curve in Fig. HJ since 
Pr depends on g*. We obtain the time lengths At/t ~ 0.015, 0.023, and 0.055. As expected, 
this approximation gives the correct value only when T gets close to T c ; otherwise, Eq. ( fl4l) 
gives just a lower bound for At. 

Fixing now an intermediate value g\ = 30 (g* pa 65), which shows more clearly the effect 
of reheating, we consider three different values of the friction (solid curves in Fig. ED, in the 
range rj ~ 0.1 — 10. These correspond to velocities which have values between v w ~ 0.1 and 




v w ~ 10~ 3 before reheating. We see that, as expected, reheating begins earlier for larger 
initial velocities. However, a variation of two orders of magnitude in fj does not change 
significantly the amounts of supercooling and reheating. For the rest of the paper, we will 
consider fj ~ 1. 
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FIG. 5: Temperature variation for gi = 30. The three solid lines correspond to v = 24QGeV and, 
from right to left, to fj = 50, fj = 5, and fj = 0.5. The dashed lines correspond to fj = 5 and, from 
right to left, to v = WOMeV and v = W~ 3 eV. 

It is interesting to examine the role of the energy scale in the dynamics of the transition. 
The model of Eq. (1231) has a single parameter with dimensions, namely, the minimum v, 
since the masses of all the particles are of the form rrii = h{V (notice that even the mass 
( I28p of the stabilizing boson is proportional to v). Thus, dimensionless quantities such as 
e.g. the ratio T c /T will not be altered if we change the value of v. This holds for all the 
quantities that are derived from the free energy (e.g., L/T 4 ,r/T 4 ), since the shape of the 
normalized effective potential in Fig. [2] is unaffected. Therefore, changing the scale v will 
not affect the dynamics of the transition, except for the expansion rate of the Universe Eq. 
( 1361) . which depends on the ratio Mp/T. 

To see the effect of such a change of scale, we have included in Fig. [5] a couple of examples 
in which the free energy is the same as before, apart from the value of v. We considered the 
QCD scale, v ~ lOOMeV, and a scale v ~ 10~ 3 eV, corresponding to a very recent phase 
transition (right and left dashed lines, respectively). Again, the temperature decreases at 
the same rate during supercooling, as explained above. However, bubble nucleation and 
reheating begin sooner. This happens because at later epochs the expansion rate H is 
slower. As a consequence, the nucleation rate T becomes ~ H 4 with a smaller amount of 
supercooling. In contrast, L/T 4 has the same value for any scale v. Thus, since 5pn is 
smaller, the temperature gets closer to T c . Hence, phase coexistence is favored in phase 
transitions occurring at later times. The parameters r and R have the same values for all 
the curves in Fig. and Eq. ([14]) yields At/t w 0.029. 

So far we have varied the parameters gi, rj, and v , which do not change the shape of the 
effective potential. We shall now consider different values of the couplings and hf, fixing 
for simplicity gj, — Qf — 10. We have checked that fixing instead and hf and considering 
different values of gj, and gf gives similar results. In what follows, we will set v = lOOMeV. 
The result is shown in Fig. El In the upper panels we plot the temperature T m reached 




during supercooling, together with the temperatures T c and T . Notice that T m is always 
closer to T c than to T . The lower panels show the time at which the temperature T m is 
reached, i.e., the duration At s of supercooling (solid line). The estimated duration At of the 
phase transition is also shown, for different values of g\. As we increase the number of light 
particles (without changing the potential), we obtain less reheating for the same amount of 
supercooling. Hence, increasing gi gives the same At s but a lower At. 
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FIG. 6: Plots of the temperatures and time intervals as functions of hb, for v = lOOMeV and 
9f = 9b = 10- Upper panels: The temperatures T c , T m , and To for hf = 0.7 (left) and hf = hb 
(right). Lower panels: The supercooling time At s corresponding to the upper panels (solid lines), 
and the total time At for different values of gi . 

The left panels of Fig. [6] illustrate the effect of a variation of hb with hf fixed. As we have 
seen in the previous section, in this case the temperature T vanishes for a value hb = Tt&i, 
where a zero-temperature barrier appears. For a value hb = hb2, the critical temperature 
also vanishes. The temperature T m lies between T and T c , so it must vanish for some value 
/i max with h b \ < h max < h b 2 (in the present case, h max m 1.3). Our numerical calculation 
does not allow us to plot the curve of T m up to this limit, because the supercooling time 
diverges for hb — > h max . This can be seen in the lower left panel. For hb > /i max the system 
never gets out of the supercooling stage. On the contrary, for hb < h max the phase transition 
completes in a finite time. Regarding phase coexistence, it occurs when At s < At. For a 
given gi, this happens up to a value of h b which is less than h max . Beyond that value, the 
supercooling temperature T m is too low for the latent heat to provide the required amount 
of reheating. Then, the estimation (|T3|) for At breaks down and the duration of the phase 



transition is just given by At ~ At s , since the phase coexistence stage is replaced by a short 
reheating (see e.g. Fig. Hj). 

If we now keep hf = hb as we increase hf, (right panels in Fig. ED, we see that T c does not 
vanish, and T decreases like l/h b as expected. In this case, the supercooling time At s does 
not diverge at any finite value of hb, and At can be considerably larger than At s . We see 
that for small values of g\ there is phase coexistence for any value of hb- On the contrary, 
for large values of gi there is no phase coexistence at all, and the estimation for At breaks 
down. The curves of At saturate for h b large because, for g t ^ 0, the parameters R and r 
cannot get close to their limits R = 1/3, r = 4/3. 

Let us now consider the case in which hf = h^, but the boson mass squared has a 
constant term pf given by Eq. (1251) so it is partially decoupled from the thermodynamics. 
The curves we obtained are similar to those in the right panels of Fig. El except that the 
temperatures meet at a finite value h m \ n . At this point the times fall to zero, since the phase 
transition becomes second-order. We have plotted in Fig. [7|the ratio At s /At for this case 
and those of Fig. El For each set of curves, the supercooling fraction increases with gi, since 
At decreases. We also see that phase coexistence is favored for hf > hb- This is because 
fermions contribute to the latent heat without enhancing the strength of the transition. 
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FIG. 7: The fraction of time of supercooling for the cases hf fixed (solid), hf = hb (dashed), and 
hf = hf, with fib (dashed-dotted). For each set of curves, from bottom to top = 35, 60 and 
120. 

In general, the energy density of radiation, p^, is much larger than the latent heat, since 
only strongly coupled particles contribute significantly to the latter. It is interesting to 
consider the case in which there are no light d.o.f. at all, i.e., gi = gu = gfi = 0. Only in this 
case the parameters can be close to the thermodynamical limits L « p^ + p\, p\ « Pi?/3. 
In the absence of light particles, all the latent heat that is released during bubble expansion 
is absorbed only by the heavy particles, which are thus more easily reheated. Consequently, 
this scenario will be the most favorable for phase coexistence. We plot the time intervals in 
Fig. El for the case hf = hb- We find, as expected, that the phase coexistence time is notably 
enhanced. For lower scales v we will have the same At but a smaller At,. 





FIG. 8: The time intervals for gi = = gn = and ht = /if,. 



V. PHASE TRANSITION DYNAMICS AND COSMOLOGY 



The cosmological implications of a phase transition depend drastically on the dynamics. 
In this section we discuss how the different steps in the evolution, namely, supercooling, 
reheating and phase coexistence, affect some of the observable products of a phase transition. 



A. Late-time phase transitions and false vacuum energy 



Late-time phase transitions have been studied in connection to the formation of large- 
scale structure and have been related, for instance, to axions, domain walls and neutrino 
masses (see e.g. [&]). In contrast to those occurring in the early Universe, which take place 
in the presence of a hot plasma with a large number of degrees of freedom, low-scale phase 
transitions happen in general in a sector with a few d.o.f. and, consequently, a small specific 
heat. Therefore, one expects a significant reheating during the phase transition, and a long 
phase coexistence stage {§]. Indeed, we have seen in section [W] that both a low v and a 
small g* favor a long phase coexistence period (see e.g. Figs H] and [5]). This stage can 
be significantly long for strongly first-order phase transitions, as shown in Fig. [6] (lower 
right panel). In particular, if all the particles have strong couplings, Fig. [8] shows that the 
coexistence of phases can last for a time At 3> tj. 

Recently, late-time phase transitions have been considered with the aim of (partially) 
solving the dark-energy problem. While the system is trapped in the metastable phase, the 
energy density of the false vacuum provides an effective cosmological constant. Thus, a false 
vacuum energy p\ ~ (10~ 3 e^) 4 could explain the observed acceleration of the Universe. 
This fact has motivated several models in which a phase transition at a scale v ~ 10~ 3 eV 
occurs in a hidden sector [Io|, U, 3l[. Such a false vacuum must persist until the present 



epoch. Hence, since the temperature of the hidden sector must be lower than that of photons, 
T 7 ~ 10~ 4 eV, the system must be in the metastable phase still at T <C v (for a discussion, 
see e.g., 0, EH). This could be achieved, in principle, in several ways, namely, due to a low 
critical temperature T c <C v lioj], due to a large amount of supercooling or due to a 
long phase coexistence stage |9( . 

The constraint on the temperature of the hidden sector, T < T 7 , comes from the Big Bang 



Nucleosynthesis (BBN) constraint on its radiation energy density, pr < 0.1p 7 . Therefore, if 
the false vacuum energy pa is to explain the observed dark energy, the temperature of the 
system must be such that 

p R {T) < 1(T 5 pa- (37) 

One possible way out of this limitation would be to assume that, although the BBN con- 
straint T < T 7 was satisfied for most of the history of the Universe, when T reached 
TV ~ 10~ 3 eV A the system entered a long phase coexistence stage at constant temperature 
[9[]. Then, as T 7 continued decreasing, the temperature of the hidden sector was stuck at 
T = T c . Thus, the BBN condition pr < p 7 would be violated only at the present epoch, 
avoiding the restriction (I3"7j) . As we have seen in section [Til a very long phase coexistence 
stage is possible. However, the effective cosmological constant during phase coexistence is 
p\ = Pa~ Pr/3 which, according to the thermodynamical bound Eq. ( fl5l) . is negative and 
does not lead to accelerated expansion. 

Due to the bound Pr(T c ) > 3p\, the condition ( 1371) cannot be fulfilled at T > T c . This 
automatically rules out any model in which false vacuum energy is dominant because the 



phase transition is yet to occur. For instance, in Ref. [10] a potential with a negative 
quadratic term — m|0 2 is considered. In that model, the temperature is assumed to be high 
enough that thermal corrections trap the system in the false vacuum. Then, the condition 
( J3"7j) is shown to be achieved for a somewhat small value of a coupling constant A. Clearly 
the thermodynamical bound is strongly violated. However, the thermal correction to the 
effective potential is assumed to be ~ T 2 2 , which corresponds to keeping only the quadratic 
term in the power expansion of the thermal integral It is then argued that the field 

is trapped at the origin as long as T 2 is large enough to cancel the negative mass squared. 
This would be correct in a second-order phase transition, in which T c = Tq. However, with 



the parameters of Ref. [10J the phase transition is strongly first-order. Hence, at T = T the 
field certainly lies in the minimum m ^ 0. The critical temperature can in fact be much 
larger than T , as shown in the upper-left panel of Fig. El where we see that T can vanish 
while T c is still of order v. 

Another possibility to attain condition (1371) is in a model with a large amount of super- 
cooling, so that Pr(T c ) ~ pa but T <C T c . This is possible in a strongly first-order phase 
transition. For the model considered in the left panel of Fig. El there is a maximum value 
of the coupling h b = /i max for which the supercooling temperature T m — > and the dura- 
tion of supercooling becomes infinite. It is not clear, however, that the required amount of 
supercooling can be achieved in a realistic model. Notice that, even when T vanishes (i.e., 
for h b = h b i), we still have T m ~ T c . In the example of Fig. El PR(T m ) w 2.37pA- Thus, it is 
necessary to go beyond h b = hj,i, i.e., to consider a model which has a barrier still at T = 0. 
In Ref. the condition (|37[) was accomplished in a specific model with T = and some 
fine tuning of the parameters. However, the thermal corrections were taken into account 
only by introducing a term ~ T 2 2 . As pointed out in Ref. [9}, this causes an unrealistically 
large value of the latent heat, which violates the constraint Eq. ( fl6l) . 



B. Electroweak baryogenesis and baryon inhomogeneities 

It is well known that the electroweak phase transition could be the framework for the 
generation of the baryon number asymmetry of the Universe (BAU). A first-order elec- 
troweak phase transition provides the three Sakharov's conditions for the generation of a 



BAU, although physics beyond the minimal Standard Model is mandatory in order to obtain 
a quantitatively satisfactory result [3;]. Due to CP violating interactions of particles with 
the bubble walls, a net baryon number density tib is generated around the walls of expand- 
ing bubbles. Assuming that CP violation is strong enough and that the baryon number 
violating sphaleron processes are suppressed in the broken symmetry phase, the resulting 
riB depends on the bubble wall velocity v w . If the velocity is too large, sphalerons will not 
have enough time to produce baryons. On the other hand, for very small velocities thermal 
equilibrium is restored and sphalerons erase any generated baryon asymmetry. As a conse- 
quence, the generated baryon number has a peak at a given wall velocity, which is generally 



KT 2 132, 33, 34 



As we have seen, reheating is always appreciable, even if there is no phase coexistence 5 . 
The temperature rise causes the wall velocity to descend significantly. Thus, baryogenesis 
is either enhanced or suppressed, depending on which side of the peak of ub(v w ) the initial 
velocity lies [H, 29]. Furthermore, baryon inhomogeneities arise due to the variation of 



Electroweak baryon inhomogeneities may survive until the QCD scale [361 . and 



affect the dynamics of the quark-hadron phase transition [36j, ]37|, [38] . The geometry of the 
inhomogeneities was studied in Refs. [1, H^]. Since baryon number is generated near the 
bubble walls, a spherical inhomogeneity with a radial profile is formed inside each expanding 
bubble. 

Notice that bubble nucleation stops as soon as reheating begins. In fact, due to the 
exponential variation of the nucleation rate with temperature, most bubbles are formed in a 
small interval 6tr around the time t m at which the minimum temperature T m is reached {23J . 
This interval is in general much shorter than the time it takes expanding bubbles to complete 
the phase transition. Therefore, it is a good approximation to assume that all bubbles are 
created at t = t m . At a later time, their walls are moving with a velocity v w (T(t)). Hence, 
all the inhomogeneities have the same profile. 

In Refs. [3. |23|. the size and amplitude of the electroweak baryon inhomogeneities were 
investigated using a simple effective potential, whose parameters were adjusted so as to 
give the desired values of the thermodynamic parameters. This approximation allows to 
vary independently parameters such as, e.g., the latent heat or the bubble- wall tension. 
These parameters, though, are generally related in a non-trivial way, which depends on the 
extension of the SM that is considered. For instance, a strongly first-order phase transition 
will have in general a considerable amount of supercooling, and also a large latent heat. 
However, the relative importance of supercooling and reheating depends significantly on the 
specific model, as can be seen, for instance, in Fig. [7J 

The amplitude of the baryon inhomogeneities, n Bmax /n Bmin , is bounded by the ratio of 
the highest and lowest wall velocities reached during bubble expansion, v max /v min . If L is 
freely varied, one can achieve values fmax/^min ~ 100 or higher [i. 23]. However, in a specific 
extension of the SM this will not be necessarily so. To examine a more realistic situation, 
we have considered extensions of the SM as in the previous sections. We find that, if we 
add a strongly coupled boson, or a boson and a fermion with hf < hb, the velocity variation 
is in general v max /v m i n ~ 1. We find a sizeable ratio only in the case in which the fermion 
dominates. The addition of strongly coupled fermions was investigated in Ref. [lj|, in order 
to make the electroweak phase transition strongly first-order. 



5 An exception could be the case of an extremely supercooled electroweak phase transition, for which 
reheating may be negligible. Such a model has been considered recently in Rcf. 35|. 



Let us consider for simplicity an extension with gf = 10 fermionic d.o.f. with mass rrif = 
hf<p, and stabilizing bosons with = gf, hb = hf, and a dispersion relation ml — ^ + hl<f) 2 , 
with fi b given by Eq. (12%]) . We obtain the plot of Fig. [HI For hf in the range of the figure 
the value of the order parameter is 4> C /T C > 1, as required by electroweak baryogenesis. The 
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FIG. 9: The ratio Umax/^min as a function of hf. 



distance scale of the inhomogeneities is given by the final size of bubbles, which depends on 
the distance between centers of nucleation. Thus, it can be roughly estimated as d ~ n^ 1 ^ 3 . 
For the present case we obtain the dashed-dotted curve in Fi g. TU1 Our results for the 
distance d agree in order of magnitude with those of Refs. 0,13^ However, we see that 
the amplitude of the inhomogeneities can be important only for small values of hf. In 
particular, fmax/^min ^ 100 is reached for values of hf for which (p c /T c < 1. Therefore, 
baryon inhomogeneities of significant amplitude are not likely produced in the electroweak 
phase transition. 




FIG. 10: The distance d between centers of nucleation for gf = (solid line), gf = gi, = 10 and 
hf = hb (dashed line), and gf = gb = 10, hf = hb, with ^ ^ (dashed-dotted). 



C. Topological defects and magnetic fields 



If a global U (1) symmetry is spontaneously broken at a first-order phase transition, the 
phase angle 9 of the Higgs field takes different and uncorrelated values inside each nucleated 
bubble. When bubbles collide, the variation of the phase from one domain to another is 
smoothed out. According to the geodesic rule, the shortest path between the two phases is 
chosen 

S3- 

When three bubbles meet, a vortex (in two spatial dimensions) or a string (in 
3d) may be trapped between them. This mechanism can be generalized to higher symmetry 
groups and other kinds of topological defects. 

If the dynamics for the phase 9 is not taken into account, the number density of defects 
depends only on the final bubble size. The probability of trapping a string at the meeting 
point of three bubbles is 1/4. Thus, the string density (length per unit volume) is ~ l/4d 2 , 
where d is the distance between bubble centers [41] . Fig. [TU1 shows the different possibilities 
for the length d. For stronger phase transitions, the bubble separation is larger, since the 
nucleation rate is more suppressed. 

Taking into account the dynamics of phase equilibration, the number density of defects 
depends also on the velocity of bubble expansion. If the latter is much less than the velocity 
of light, the equilibration between the phases of two bubbles may complete before a third 
bubble meets them, thus reducing the chances of trapping a string. Consequently, reheating 
hinders the formation of topological defects. 

In the case of a gauge theory, a spatial variation of the phase 9 is linked to a variation of 
the gauge field jUJ . As a consequence, a magnetic field is generated together with the phase 
difference in the collision of two bubbles. Then, one can say that a vortex is formed whenever 
a quantum of magnetic flux is trapped in the unbroken-symmetry region between three 
bubbles. The phase equilibration process is thus related to flux spreading, and depends on 
the conductivity of the plasma. Bubble collision constitutes also a mechanism for generating 
the cosmic magnetic fields (see e.g. [42j). This mechanism may take place at the electroweak 
phase transition, where unstable cosmic strings and hypermagnetic fields may be formed. 
The latter are subsequently converted to U (l) em magnetic fields. 

A detailed calculation of the density of defects and the magnitude of the magnetic fields 
is beyond the scope of this paper and we leave it for future research. Although some 
simulations have been made (see, e.g., j43||). several simplifications are generally used, which 
include assuming a constant nucleation rate and a constant bubble wall velocity. As we have 
seen, this situation is hardly realistic. Moreover, the formation of topological defects and 
magnetic fields depend strongly on the dynamics of the phase transition. In particular, a 
long phase coexistence stage with a very slow bubble expansion will affect significantly the 
mechanism of phase equilibration during bubble percolation. 



VI. CONCLUSIONS 



In this article we have investigated the different stages in the development of first-order 
phase transitions of the Universe. In particular, we have studied the amounts of supercooling 
and reheating. If the entropy discontinuity As(T c ) is larger than the entropy decrease 
5 s = s(T c ) — s(Tjv) during supercooling, a phase-coexistence stage is reached. Then, the 
total duration of the phase transition can be calculated analytically. The ratio At/(2H)^ 1 
depends only on the parameters r = L/pn and R = pa/pr. If As(T c ) < 5s, supercooling 
lasts for a time which is longer than At. In this case, there is no phase coexistence, and 



At gives only a lower bound for the total duration of the phase transition. We have shown 
that thermodynamics constrain these parameters to the region R < 1/3, r < R + 1. These 
constraints should be taken into account when the dynamics of a particular phase transition 
is considered, since approximations for the effective potential may violate them, and thus 
the analysis may lead to incorrect results. 

With the help of a simple model, we have analyzed numerically the role of different 
parameters in the dynamics of the phase transition. We have verified that phase coexistence 
is more likely in later phase transitions, since both a lower energy scale and a smaller number 
of degrees of freedom favor reheating. In addition, we have seen that changing the viscosity 
of the surrounding medium does not affect significantly the dynamics of supercooling and 
reheating, although it affects the velocity of bubble walls. The incorporation of bosons to a 
given model strengthens the phase transition, so the effect on the dynamics is to enlarge the 
latent heat and suppress the nucleation rate. As we have seen, the latter effect is in general 
stronger, so adding bosons favors supercooling. On the contrary, adding fermions in general 
weakens the phase transition and at the same time increases the number of d.o.f. We have 
checked that in this case phase coexistence is favored. 

We have studied how our general results on phase transition dynamics may affect some 
of the cosmological consequences. For instance, in the case of dark energy from a phase 
transition, we have shown that the thermodynamical bounds rule out some models. Besides, 
we have analyzed the effect of dynamics on two important parameters, namely, the number 
density of bubbles and the amplitude of the velocity variation during reheating. As we 
have seen, these quantities are relevant for the generation of different cosmological relics, 
e.g., baryon inhomogeneities, topological defects and magnetic fields. In particular, we have 
found that it is difficult to obtain baryon inhomogeneities of sizeable amplitude in realistic 
models of the electroweak phase transition. 

We believe that our results on the dynamics can be applied to a wide class of phase tran- 
sitions of the Universe, and the discussion on the cosmological consequences can be extended 
to several interesting possibilities, such as, e.g., the formation of baryon inhomogeneities in 
the quark-hadron phase transition j^] or the generation of gravitational waves (6j. 
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APPENDIX A: APPROXIMATIONS FOR THE THERMAL INTEGRALS 

In this appendix we consider expansions of the functions I±{x) for small x and large x. 
The integrals in Eq. (j22l can be evaluated numerically. However, a numerical computation 
in the effective potential increases significantly the total computation time. Indeed, notice 
that for each temperature, we must find the minimum (f> m (T) to compute several quantities 
derived from J-"(T). Moreover, the calculation of the bounce action Sa(T) requires the time- 
demanding overshoot-undershoot technique to solve Eq. ( 13T1) for the bubble profile at each 
T. Therefore, it is useful to employ analytical approximations for the thermal integrals. 



Following the derivation of Ref . 44| , we can obtain the expansions of I± (x) in powers of 
x. For bosons we have 
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where a fe is given by log a;, = 3/2 — 27 + 21og(47r) , with 7 the Euler constant; ( is the 
Riemann zeta function, and T is the Gamma function. The expansion for fermions is 
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where aj is given by log a/ = 3/2 — 27 + 21og7r. For any value of x we can get the desired 
precision by keeping enough terms in these expansions. For example, keeping up to I — 5 in 
I and I = 12 in I + , we obtain a precision of 10~ 8 for < x < 2. 

The expansion for large x can be obtained by changing the variable of integration to 
z = \Jy 2 + x 2 and expanding the logarithm in Eq. (1221) in powers of e~ z (see Ref. 18]). 



k=l 



I T ( x ) = - ^ ^ / dz zVz 2 - x 2 e' kz . (A3) 



For each k, the integral yields x 2 Ki (kx) /k, where K2 is the n = 2 modified Bessel function 



of the second kind 45] K n (z). Hence, we obtain the expansions 
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Notice that the integrals in Eq. ( 1A3I) are of the order of e~ kx , so the terms in this expansion 
decrease with powers of e~ x . Therefore, in general we will obtain the desired precision by 
considering a few terms. For example, for x > 10 we obtain AI/I < 10~ 10 by keeping only 
the first two terms in (1A4|) . For x > 2, keeping terms up to k = 7 in the expansion gives a 
precision AI/I < 10~ 8 . As a rough estimation of the error of the truncated expansion, we 
note that the fc-th term is ~ x 2 e~ kx /k 2 , and the error is given by the ratio of the (k + l)-th 
term to the first term, AI/I ~ e~ kx /k 2 . 
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